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I. INTRODUCTION 


Basic to the design of reflecting antenna systetns Is the need to 
compute the far-fleld radiation pattern for a given antenna geometry, feed 
placement, and feed radiation characteristic. A simple, yet fast solution 
to this classic problem has proven to be a formidable challenge. This Is 
especially true for very large reflectors with extensive current distri- 
butions and off-focus feeds with their associated asymmetry. These problems 
are made troublesome because of their large computational requirements. 

Recent Interest In this problem [1-4] has resulted In the new appli- 
cation of powerful mathematical techniques to produce a more time-efficient 
means of computing the radiation Intensity. Galindo- Israel and Nittra [l] 
have reformulated the two-dimensional radiation Integral to a double Fourier 
integral “representation." The Fourier Integral Is evaluated by expanding 
the Integrand In a double sunmatlon of circular functions and modified 
Jacobi polyncmlals. The result Is a rapidly convergent series representa- 
tion of the radiation Intensity, the zeroth term being Ji(x)/x, the 
pattern for a uiilform, cophase aperture. Mittra, et al. [2] extend the 
work of the first paper to offset paraboloid reflectors and Introduce a 
more rigorous corrective technique to compensate for the .lOn-Fourler 
structure of the radiation Integral. 

This paper presents a new method which evaluates the radiation 
Integral on the antenna surface. The algorithm Is computationally efficient 
and produces the far-fleld radiation pattern along planer cuts at any angle 
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I. INTRODUCTION 


Basic to the design of reflecting antenna systems Is the need to 

* 

compute the far-f leld radiation pattern for a given antenna geometry, feed 
placement, and feed radiation characteristic. A simple, yet fast solution 

to this classic problem has proven to be a formidable challenge. This Is 

* 

especially true for very large reflectors with extensive current distri- 
butions and off-focus feeds with their associated asymmetry. These problems 
are made troublesome because of their large computational requirements. 

Recent Interest In this problem [1-4] has resulted In the new appli- 
cation of powerful mathematical techniques to produce a more time-efficient 
means of computing the radiation Intensity. Galindo-Israel and Hittra [l] 
have reformulated the two-dimensional radiation integral to a double Fourier 
integral “representation." The Fourier integral is evaluated by expanding 
the integrand in a double summation of circular functions and modified 
Jacobi polynomials. The result is a rapidly convergent series representa- 
tion of the radiation Intensity, the zeroth term being J^CxI/x, the 
pattern for a uniform, cophase aperture. Mittra, et al . [2] extend the 
work of the first paper to offset paraboloid reflectors and introduce a 
more rigorous corrective technique to compensate for the non-Fourier 
structure of the radiation integral. 

This paper presents a new method which evaluates the radiation 
integral on the antenna surface. The algorithm is computationally efficient 
and produces the far-field radiation pattern along planer cuts at any angle 
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f through the thrce-d linens 1 one 1 pattern. The physical optics approxima- 
tion Is used to compute the, Induced surface current which Is the input 
to the algor Ittm, The method Is developed for focal -point and translated 
feeds and Is easily extended to offset antennas and arbitrary surfaces. 

On the premise that simplification and efficiency may come from 
first generalizing* the radiation Integral Is reformulated to three dimen- 
sions. The result Is shown to be a triple Fourier Integral. To evaluate 
this Integral, an extremely fast algorithm Is Introduced which evaluates, 
along planes of constant i>* a subset of the total three-dimensional Fourier 
transform results. No approximations are made other than those normally 
associated with digitization and the discrete Fourier transform (DFT). 

To further reduce the computation time, the Winograd Fourier transform 
algorithm (WFT) is used In place of the standard radix-2 FFT when a OFT 
is called for in the algorithm. The any ^ angle feature of the program 
is implimented using a technique similar to one used In computerized 
tomography (cross-sectional ,x-rays) [5]. 

While other methods exist for evaluating the radiation Integral, 
this new theory brings a fast, simple and direct approach. Due to its 
speed it Is especially useful for very large asymmetric antennas. 

II. THE THREE-DIMENSIONAL RADIATION INTEGRAL 
AS A FOURIER TRANSFORM 

The far-field radiation intensity of a volume distribution of 
current may be expressed in terms of the three-dimensional radiation Integral 




( 1 ) 
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The geometry for (1) Is given In Figure U In formulating this equation# 
It Is assumed thst the phase center Is at the origin# the observation 
point Is far from the current distribution and the currents are bounded 
In a region dimensionally small compared to R. 

The -r'*R (R a unit vector) term Is the distance from the source 
point to a plane through the origin and normal to R. Hence# the exponent 
k(r '-?'•&) Is the total phase delay modulo 2itR. 



Figure I. Geometry for the radiation Integral 
Using the direction cosines u, v# and w. we can write 
r'»R ■ r^w + r^u + r^v 

■ r'cose'cose + r'cosf*sine’cos<*s1n0 + r'sln^'slnd'sln^slne 
« r'cose'cose + r*slne‘s1necos(<>-^*) 


The radiation Integral of (1) may now be rewritten as 


4*0 


( 2 ) 


^^jkr'(s1no*s1n'4»cos(M’)Vi.asine'dr'do’f:4' 

Equation (2) 1s a three-dimensional Fourier transform In spherical coordi- 
nates as presented by Bracewell [6]. An Important corollary to this 
observation Is that the far-fleld radiation pattern may be computed as 
the three-dimensional Fourier transform of a volume current distribution. 
Techniques will be presented for efficiently evaluating equation (2). 

If the current exists only on a surface within the volumet the 
volume current function may be expressed In terms of a surface current as 
K(r‘,0',4') * 7(t/*,4')fi(r’-p) 

where p * p(e*,^') defines the geometry of the surface. Substituting this 
(j{i;Kpress1on into (2) and integrating with respect to r' (with the aid of 
the sifting property of the Dirac function) we arrive at 

E(e.*) . /” r 

♦•0 e»o ( 3 j 

^ejkp[s1ne's<necos(*-*')]p2sin9'de'd4.' 

Equation (3) agrees with equation (5) of Galindo- Israel [1], To 
phrase it In exactly the same form requires two modifications. First, 
the phase center must be moved from the origin to a point defined by the 
feed position vector t. Secondly, the integration must be transferred 
from the reflector surface to the aperture plane. The latter Is accom- 
plished by use of the "Jacobian of the transformation" between the reflector 
surface and the projected aperture. A single component of J(e',4'), say 
«lj((o',4')t may be transformed to an "equivalent" aperture distribution 
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over the circular disk of the aperture. This results In the transformation 

+ (dz/ds)?'sdsdi‘ 

« f(s.^*)sdsd^' 

f(s,^') is the equivalent aperture distribution and Is not to be confused 
with the physical distribution on the aperture such as might be derived 
approximately by the ray technique. 

Equation (3) now becomes 

E(e.*) . /** t 

♦•0 $■» (4) 

,jjkB[s1ne's1necos(*-*')]jj5j^. 

Which is the same as (5) in 6al1ndo*israel [1] except for a difference 
In phase center. Equation (4) Is not a Fourier transform as 
exp(jkp cose 'cose) Is a function of both source and observation coordi- 
nates. 

m. AN EFFICIENT ALGORITHM FOR THE THREE-DIMENSIONAL 
DISCRETE FOURIER TRANSFORM 

The three-dimensional Fourier transform of a volume source Is 
defined by Bracewell [6] as 

F(u,v,w) « r r r f(x,y, 2 )e’’^^^^“^ + vy + (5) 

•00 -ee ««D 

where u» v, and w are the direction cosines 
u • slnecosi^ 

V ■ sinesin^ (6) 


w > cose 
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Note that while there appear to be threi? Independent variables in this 
transform* i.e., the frequencies u, v, and w, there are in fact on«y two 
geometric variables, e and Ihese may be solved for as 

e • sin"Vu^ t V* ■ cos"^w 

(7) 

a ■ tan“*v/u 


Similarly, the three-dimensional discrete Fourier transform (DPT) 
is given by the triple sunmatlon 

2 

( 8 ) 


Ns"! Nj-1 N3-I j2it_ ^ j2w_ t. j 2if^ ^ 

FCkj.k-.k,) . I I 

n2*0 ni*0 n3*0 


for kj * 0. 1, 2, ... , 

* 0, 1, 2, . . .. , Ng”! 

^3 * 8 , l| 2 f . • . , Ng-l 


This procedure may be broken into three operations by the use of inter- 
mediate results. Consider the intermediate transforms 


N 3 -I ' I, 

g(k3,n9,n,) « I f(n3,n2,n,)e“^R3 

* ns-O 

nj *< 0, 1, 2, . . . , Nj-1 

n^ * 0, 1,2, . . . , N^— 1 

k 3 * Of 1 , 2 , ...» N 3— 1 

Ni-1 j2tt 

h(k 3 .nj,.kj) « I g(k3.n2.n,)e"^R7”i'‘i 

ni*0 

kj ^ 0, 1, 2, . . . , Nj-1 

* 0, If 2, . . . , Ng— 1 

. W 3-1 


(9) 


( 10 ) 


k3 « 0, 1. 2, 


• • • 
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From these the DFT is found to be 

FCkj.kj.k,) - Ii(k3,n2.kiit’'^^*''’ 

kj • 0. 1, k Nj-1 

kg • Of It 2* • • » t Ng"! 

kg ■ Of 1» 2.p » • » t ^3"^ 


( 11 ) 


Comparing the frequency viriables of the Fourier Integral and the 
DFTt we note the correspondence 



where each T Is the sampling period for the respective axes. Combining (12) 
and (7), we arrive at expressions for 9 and 4 In terifis of the frequency 
Indices of the DFT 


4 • tao 


A direct evaluation of the DFT 1n ( 8 ) Is impractical, even for 
modest transform lengths. Nor does use of an FFT algorithm provide suffi- 
cient Improvement. However, the computation time may be reduced by orders 
of magnitude in the ca«;e where one is Interested in a limited amount of 
the total DFT information. It is possible in this way to compute a two- 
dimensional planer cut through the far-field radiation pattern with no 
compromise in accuracy. Once the algorithm is developed, it may be extended 
to Include an any angle planer cut. 

The simplification Is begun by selecting ■ 0. From equation (13) 
this requires that 4 ■ 90^. Next, the value of kg Is chosen. From (13) 



e jj 

thli determines the poler anyte • • iin‘*(k2/NtT2) end forces 1(3 to • 
specific value piven by 

Clearly* there are going to be digitization problems in evaluating ( 14 ) for 
K3 based on a selected Integer value of k^. This works to our advantage 
in decreasing the computation time. 



Figure 2 . Antenna Geometry for the DFT 


The summation in equation ( 9 ) may now be reduced. k3 is a deter- 
mined single value and along the 03 dimension (see Figure 2 ) there Is 
only one non-zero data element for each (n^.n^) pair. Hence ( 9 ) becomes 


g(k3,n2»nj) 


f(n3,n2.ni)e 


-% 3 k 3 


( 15 ) 


n^ * Of If 2* ... • Hj-l 
n^ * Of lf2f . • • f 

k3 a constant 
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Equation (IS) describes a coalescing of the data to a planer array and 
hence a reduction to two-dimensions as shown In Figure 3. Each point In 
the 0^02 plane corresponds to a g(k3»n2ini). The fCnsinaini)* (I.e. the 
current density) data will be present In polar form so that the complex 

9 . 

multiplications Indicated here are simple additions of to the pha$e 

9_ 3 

of f(n3,n2.ni). Furthermore, the phase term ^3^3 need be calculated 
only once for each 03 and recalled as a linear array variable for the 
Individual phase additions. 



Figure 3. The coalesced two-dimensional array 
With the selection of ■ 0 comes the simplification of equa- 
tion (10) to 

Ni-1 

h(k3.n2,0) * I g(k3,n2fni) (16) 

nj«0 

ng ^ 0, 1, 2, ... I N 2 -I 

k3 a constant 

Equation (16) describes the coalescing of a two-dimensional array to one 
dimension. The g(k3,n2»ni) data along each row are summed to form the 
linear array hCka.nj.O) In n2. 

The final operation is described by equation (11) 
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FCka.kj.O)- Ii(ks,n2,0)t‘^^2''2 (17) 

n2*0 

kj ■ NjTj^-ikj/NjTj)^' 

Note that (17) Is not a DFT since the sum Is performed for ouly one value 
of kg. It Is In fact only one of the Ng harmonic constituents of the 
OFT corresponding to one value of e in the spectrum. In practicet one 
may want to compute the DFT of h(k3,n2*0), as digitization problems In 
evaluating kg make (17) applicable for many values of 6 (or k 2 ) based on 
a single value of kg. 

Generally, it Is not necessary to execute the above process for 
more than several values of kg. For NT * (2520) (1/8) In equation (14), 
kg ■ 315 for kg In the range of zero to 17. This requires that for e 
between zero and 3.1^, the generation of the h(k 3 ,n 2 > 0 ) one-dimensional 
array by coalescing to a plane and thence to a line need be done only once. 
Radiation data at e Increments within this bound are found by evaluating 
(17) for kg between zero md 17. Thif* may be done by computing the FFT 
or WFT [7] of h(k3,n2.0) and using only the first 18 transform results. 
Subsequent integer values of kg yield e bounds of 3.1^ to 5.5^, 5.5^ to 
7.1®, etc. 

The three-dimensional DFT produces a spectrum In the three variables 
kj,kg, and kg. Such a thing is difficult to represent graphically and 
fortunately, in this case it is not necessary as k^ « 0. What results, 
then, is a two-dimensional spectrum or a surface in kg and kg. Only a 
portion of this surface is of interest; that section defined by the elliptic 
relationship between k 2 end kg In equation (14). It is the spectrum 


n 

•long this elliptic arc (see Figure 4) that describes the far^fleld 

(• 

radiation pattern for a constant f cut. 



Figure 4. Two-dimensional spectrum for ^ 

The DFT described In (8) Is a "one-sided" transforma t1on» i.e., 

It operates on data sequences in positive time or position only. Further- 
more, the DFT is cyclic and requires that the input data be periodic. 

The first condition leads to a necessary format of data Input. The 
second requires a large number of zeros to be embedded in the data to 
adequately isolate the antenna from distant images regularly spaced three- 
dimensionally about it. The geometric association of data to the (ni,n 2 »n 3 ) 
Indices is demonstrated in Figure 5. The location of the origin for the 
n^ anJ axes is along the center line of the antenna and is a clear ex- 
tension from two-dimensional DFT theory. There is no obvious point of 
s^imetry along the n 3 axis to fix the three-dimensional origin. Conse- 
quently, it would be nice to discover that the choice did not affect the 
DFT results and could be made arbitrarily. It will be shown that this Is 
true for the amplitude results, but not the phase. 
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Figure 5. Data Input format for the three-dlmensilonal DFT 
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Equation (IS) contalnr the coordinate dependence of the DFT. A 
shift In the origin along the n 3 axis by an Integer amount y will alter 
the phase term to 

When this data Is coalesced to a line, the resulting one-dimensional data 
vector will have a constant phase shift as seen from equation (16) 

nj«0 

■ g(k3tn2.ni) 

ni*0 

« e"’^^h(k3,n2,ki) 

Hence, a constant phase shift. 6. will be present In all the DFT terms, 
but the amplitude results remain unaltered. A constant phase shift Is 
certainly not a problem. Hovvever. difficulty does arise when the algorithm 
Is repeated for other values of k 3 . Recall that repeated application is 
necessary in order to sweep through a broad range in e and that the number 
of separate values of is determined by this e range (three applications 
for e out to 7.1°). Clearly, e is a function of kg and the separate com- 
posite ranges will each be shifted by a different constant phase. If phase 
Information is important, it is an easy matter to compute these phase terms 
and subtract them from the results. 

Finally, the direction of the ng coordinate is Important. The 
Fourier transform may be formulated with either a plus or minus exponential 
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phase tehn in the integrand. The n^ coordinate orientation shown in 
Figures 2 and 5 is consistent with the negative phase form and when used 
with this definition produces correct transform results. 

IV. AN ANY ANGLE PATTERN CUT 

The selection of kj ■ 0 in the previous section was significant 
in that it allowed us to simplify the three-dimensional DFT by coalescing 
the planer data to a linear array. However, it also limited our considera- 
tion to the * 90° plane. It is imperative that we be able to examine 
the radiation intensity along other planer cuts through the pattern. This 
may be done using a projection technique similar to that used in computer- 
ized tomography [5J and multiangular scanning in gases [8]. 

Equation (16) indicates that the planer data (see Figure 3) be 
summed along rows of constant to produce a one-dimensional data vector 
along the n£ (♦■90°) axis. Following the same procedure, the data may be 
projected to a diagonal line n^ (see Figure 6) at an angle ♦ to the n^ axis. 
Conceptually, this is equivalent to a rotation of the coordinate system 
about the 03 axis. The resulting data vector is then operated on by 
equation (17) to produce a planer cut at the new angle 

From Figure 6, the complex data arc summed along grid tubes (of 
width t) perpendicular to the diagonal axis. Looking along these tubes, 
one observes that the date points do not lie in the centers of the new 
grid squares and that a few squares contain no date points, while others 
contain two. The first difficulty is related to quantitization error 
which is always present when converting a continuous function to discrete 
data. The assumption that an interval of a continuous function may be 




Figure 6. Coalescing to a diagonal line for an any angle cut 
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represented by a constant value Is a first order type of approximation; 
that the constant is not the value of the function at the center point is 
second order. 

As for the second difficulty, the concern here is that the sum along 
each grid tube be essentially the same as if resulting from a regular grid 
with one and only one value associated with each square. This will be true 
if the sampling rate Is sufficient and If the sums are adjusted to account 
for any Irregularity In the number of data In each tube. 

V. CALCULATION OF THE INDUCED ANTENNA CURRENTS 


The physical optics approximation was used to compute the antenna 
surface currents. This comes from applying the magnetic boundary condi* 
tions 

J « 2n K IT (18) 

Generally, this Is considered sufficiently accurate for studies of the 
main beam and several side lobes. It is only necessary to determine the 
unit vector (n) normal to the surface and the magnetic field (it). The 
currents are derived for a paraboloid reflector with a (-cose')” offset 
feed and a feed displacement = (x^,yg, 2 ^). Assuming a unit electric 
polarization in the y-di recti on, a l/p* space divergence of the field and a 
phase delay, we determine 
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The geometric variables for (19) are defined In Figure 7, 



/ 

Figure 7. Geometry of focal -point and translated feeds. 

VI. COMPUTATIONAL RESULTS 

Performance of the new three-dimensional algorithm Is verified by 
comparing computed results with other published work. A driver program 
was used to generate, on the paraboloid, a surface current of constant 
amplitude. Furthermore, a phase advance assigned to the 

current at each point so that the field at the aperture would be both uni- 
form and cophase. This provided a test of the new algorithm with classical 
theory. The computed result was the expected J|(x)/x (see Figure 8). 

Radiation patterns were computed for currents derived by the physical 
optics approximation. The reflector is characterized by an f/D » 0.5, and 
the feed pattern as a circularly symmetric (-cose')^. The value of n was 
chosen to produce a -lOdB taper to the edge of the reflector. Figures 9(a) 



Polar Anglo (dogroeo) 

Figure 8. Radiation pattern for a uniform, cophase aperture 
by the three-dimensional DFT 
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tnd 10(a) present radiation patterns In the 90® plane for a focal - 
polpt feed and an Offset feed (y, « -l.25x). The phase results are given 
In Figures 9(b) and 10(b). These results may be compared with Figures 

3 and 4 of Galindo- Israel and Mittra [1] t^here amplitude and phase patterns 
were computed for similar antenna parameters u$1ng the Jacobi polynomial 
method. Careful examination reveals strong agreement between the two methods. 

A more rigorous test of the algorithm Is to translate the feed 
1.25X along a line bisecting the -x and -y axes and compute a 4 ■ 45® cut 
through the pattern. In this case, the feed displacement parameters are 

* *0.8838831 and * 0. The resulting pattern should be essentially 
the same (for e<10®) as computed for y, ■ -1.25x, x^ ■ » 0, and ^ ■ 90®. 

Excellent agreement is clearly seen In Figure II. Due to polarization, the 

a ■ 45® results fall Increasingly further below the 4 ■ 90® results as e 
becomes larger. 

Figure 12 demonstrates the flexibility of this new algorithm. For 

► 

a feed offset y. ■ -1.25x, and x^ ■ z^ « 0. pattern cuts are computed for 

6 C C 

4 ■ 90®. 60®. 30®. and 0? Repeated application to other 4 angles can 
adequately describe the three-dimensional nature of the radiation pattern. 

By experience, the new three-dimensional algorithm was found to be 
computationally very fast. Timing was done on the Goddard Space Flight 
Center's IBM 360/91 using the interval timer available in the system 
library. With a very close x/8 samplingr the epu run time per pattern 
(for the general case) was approximately 25 seconds. A x/4 sampling was 
found to be sufficient to study the main beam and several side lobes (see 
Figure 13), and ran In a mere seven seconds. Circularly symmetric current 
sheets may be treated as a special case and consequently, run times re- 
duced by about 75 percent. For the epu times cited, most of the time Is 




POLAR ANGLE (DEGREES) 


Figure 9. Gain and phase patterns for a focal-point feed. 
D/X • 50, f/0 ■ 0.5, t • x/8 and N ^ 2520. ' 
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T - X/Z and N « 2520. 




ANTENNA GAIN (OB) 
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POLAR ANGLE (DEGREES) 

Figure 13. Comparison of X/8 and a/ 4 sampling for an offset feed problem. 
O— OT * a/8; a T « a/4. 
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consumed In computing the Induced currents. The pertlal three-dimensional 
DFT results are computed In much less than half of these times. 

Vn. CQtICLUSIONS 

A method of evaluating the radiation Integral on the curved surface 

% 

of a reflecting antenna has been presented. The result Is a two-dimensional 
radiation cross-section along a planer cut at any angle ^ through the far- 
fleld pattern. This section is produced by evaluating the radiation 
integral via a three-dimensional Fourier transform, A unique feature of 
this method is a new algorithm for evaluating a subset of the total three- 
dimensional DFT results. The algorithm is extremely fast so that the 
computer time required to produce a radiation pattern Is primarily determined 
by the computation of the antenna currents. 

High quality gain and phase results have been computed for a parabo- 
loid reflector with translate feed. However, the method Is easily extended 
to offset antenna systems and reflectors of arbitrary shapes. The new 
method provides a direct but fast approach to the analysis of large asymme- 
tric reflector antennas. 
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